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The agent-based model of stock price dynamics on a directed evolving complex net- 
work is suggested and studied by direct simulation. The stationary regime is maintained 
as a result of the balance between the extremal dynamics, adaptivity of strategic vari- 
ables and reconnection rules. The inherent structure of node agent "brain" is modeled 
by a recursive neural network with local and global inputs and feedback connections. For 
specific parametric combination the complex network displays small-world phenomenon 
combined with scale- free behavior. The identification of a local leader (network hub, agent 
whose strategies are frequently adapted by its neighbors) is carried out by repeated ran- 
dom walk process through network. The simulations show empirically relevant dynamics 
of price returns and volatility clustering. The additional emerging aspects of stylized 
market statistics are Zipfian distributions of fitness. 
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1. Introduction 

The design of the decentralized systems of autonomous interacting agents with 
abilities to automatically devise societies so as to form modules that accommodate 
their behavior via social and economic norms in emergent ways is highly challenging. 
Agents in such groups plan and act in a way to increase their own utility. 

As a general theoretical framework for starting of the statistical consideration 
of interacting agent entities we take into account the Ising model. Being made 
simply by binary spin variables, the model is able to reproduce different compl ex 
phenomena in different areas of science like biology sociology ^ economy 121415] or 
informatics 

Looked at from the perspective of the economics this example has a great impor- 
tance because it demonstrates that a basic interaction between the spins (agents) 
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can bring non-trivial collective phenomena. The parallels between fluctuations in 
the economic an d mag netic systems afford an appli cation of spin models to the 
market statistics I 7 | 8 | 9 | -j>he attempts | 10 | H | 12 | 13 | 14 | examme a context with the 
minority game theory Basic terms borrowed from the magnetic systems have 
been built in to the spin agent concept: exchange interaction between spins and 
interaction with random field. The approaches based on the minority game assume 
competition of the short-range ferromagnetic and global magnetization terms that 
are crucial for modeling of bubbles and crashes. The justification for our present 
formulation comes from^^. It assumes that each among i = 1, 2, . . . , L interacting 
traders (agents) owns her/his regular lattice site position i and a corresponding spin 
variable S^\i) G { — 1, 1}, where the upper index (t) labels the market time. Each 
agent, has an attitude to place buy order S^(i) — 1 or to place sell order described 
by S^(i) = — 1. The variable S^(i) is updated by an asynchronous stochastic 
heat-bath dynamics expressed in terms of the local effective field. 

The collective effect of spins is characterized by the instant magnetization 



that is interpreted as measure price imbalance. Therefore the logarithmic price 
return can be written as 



where A is the liquidity constant. The predominance of buy orders manifests it- 
self throughout > 0. In that case the instant price p( t+1 ) exceeds pw. In an 
analogous manner nv-*' < describes the fall of the stock price. 

In the spin models the non-Gaussian contribution to distribution of magnetiza- 
tion is formed due to spin-spin interactions. The realistic models of social/economic 
networks can serve as topological substrate for spreading of the interaction effects. 
Our recent approach to spin market models has been formulated for network ge- 
ometry. The attempt to postulate dynamics of network that is coupled to the spin 
degrees of freedom is analogous to structurally dynamic cellular automata in 
which the conventional cellular-automata rules are generalized to formalism where 
geometry and matter are dynamically coupled. 

Since our former attempt did not give satisfactory topological results, one of the 
aims of our present work is to propose richer and more reliable formulation. 

The main interests of our present research is focused on: 

(1) An alternative protocol of the network reconnection assuming the network 
segmented by local " authorities" (local leaders - hubs) . 

(2) The suggestion of rules that support formation of the "small world" and 
"scale-free" paradigms for social networks. 

(3) The including of the small-scale intra-agent cognitive-like sensorial struc- 
tures sustained by interactions within the complex network of agents. 




(1) 



In p^+D/pW 



= m 



(2) 
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The plan of the paper is as it follows. In the next section [2] we discuss the basic 
network properties. The overall dynamics are given in section [3l A more detailed 
definition of the model items appears in subsections 13 . Hl3~5l In section [4] we present 
statistical characteristics extracted from our simulation. 



2. The network topology 

Let's suppose that the market structure is determined by the underlying complex 
network defined by dynamical rules for active links between agents. Consider the 
directed network (graph) of labeled nodes T = {1, 2, . . . , L}, where node i G T at- 
taches via N out directed links to its neighbors X n (i) G T, n G 7 out ={1,2,..., -/V ou t}, 
i.e. the graph is -/V out -regular. Two outgoing links Xi(i) = 1 + imodL, ^(i) = 
1 + (L + i — 2)modL of node i e I create the bidirectional cycle static subgraph (L- 
gon). The reconnection rules are applied exceptionally to the links Xn\i), n G I ut 
= {«i;3 < ri\ < N out ] n% G /out}- The using of static module ({1,2} C I ut) 
guarantees the preservation of network connectedness at any stage t. 



3. The formalism of co-evolutionary dynamics 

Formally, the stochastic co-evolutionary dynamics of agents can be described by the 
recursive formula 

n (m) =u ( n (t) ) (3) 

including the composed configuration 

TT (t) = {n (t )(i),n( t '(2),...,n< t )(L)} (4) 

that consists of single-agent particulars 

intranet intra-agent spins IT2 (i) = {s^(i, fc)}feer int r 

n (t) (i) = I strate S ic variables n j (0 = Uintri^^l)}, ( 5 ) 

where k, q G Fi ntr 

network links Tl^ (i) = {Xn\i)} , Tl G 7 ut • 

The intranet of agent i includes fully connected recurrent network consisting of 
nodes Fi ntr = {1, 2, . . . , 7Vi ntr } nodes [for more details sec modification Eq. fTTj) ]. 
The nonlinear operator U entails the overall effect of the following single-agent 
operators 

local field U ss (i) acting on Llgs'' (i) 

&Ad(«a) acting on IIj (i a ) ^ 
reconnection t^R c (i r ) acting on IIy (*r) 

extremal dynamics [/Ex(«minF) acting on IlH \i m inF) ,ny^(i m i„F) 
where i m i n F is defined by minimum of fitness [see Eq. (|13p in further] 

F(w F ) = minF(i). (7) 
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The pseudo-code corresponding to dynamics Eq. ^ is described in the next five 
steps: 

• loop I over the Monte Carlo steps per network node 

loop II over the network of agents 

1. pick agent i randomly 

2. perform the field and spin update U ss (i) 
[see Eqs.® and © of subsection |3~T] . 

3. pick two agents randomly: 

i a for the adaptivity rule and 
i T for the reconnection rule. 

4. apply C/Ad(«a) 

with probability P\d (see subsection 13. 3p 

5. apply U Rc {i r ) 

with probability P Re (see subsection 13. 4p 

end of loop II 

apply the extremal dynamics J7ex(Wiif) to agent i m i n F having 
the smallest fitness within the entire system, 
end of loop I 

After it the following steps are carried out: 

• (i) store instant m calculated according Eq.fTT]), 

• (ii) calculate the price return Eq.([2|) and volatility as \m^\ 

• (iii) update fitness F(i) for all agents according Eq. fTB")) (see subsec- 
tion^ • 

3.1. U ss : decision making via local field, interface between extranet 
and intranet 

Through links X n (i) agent i gains the game-relevant information about the external 
nearest world. For directed network topology we define the integral spin value 

Snn(i)^^- £ S(X n (i)). (8) 

out nei out 

The intra-agent idiosyncratic structure is modeled in a abstract way that resem- 
bles the modular organization of the routing control unit The present variant of 
the model goes essentially beyond the elementary single-spin intra-agent descrip- 
tion. We have looked in this direction by considering trader's states encoded by 
Ising spins (neurons) coupled by the weighted links of the fully interconnected in- 
tranet. We suppose the agent's i architecture consisting of -/V; n t r spins {s(i, &)}/cer intl j 
s(i,k) G { — 1,1}. In further, we distinguish between the inner s(-, •) and extranet 
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S(-) symbols. The effective local field h\ oc (i, k) is suggested in the form 

h\ oc (i, k) <— h stoch (i, k) + — — V" J hl t r (i,k,q)s(i,q) , (9) 

nntr t 

q=i k 

q- ke r intr 

where /i s toch(*j k) is the Gaussian stochastic variable (i.e. the source of local muta- 
tions) and Ji„t r (i, fc, g) is the M ntr x -/V; ntr system of the intra-agent pair couplings 
(weights); the term l/(-/Vi ntr — 1) is included for normalization reasons. Three inputs 
of intranet are supposed: 

s(i,3)<-m, s(i,4) <- 5(i) , s(i, 5) <- 5 nn («) ■ (10) 

Two of the spins s(i, k), k = 3, 5 are linked to outer complex network represented 
by m and 5' nn (i). The feedback to s(i,4) is related to self-control. The state of fcth 
unit of agent i is recalculated according 

s(i,k) «- sign(/iio C (i,A;)) , when fc € r intr = {1, 2} U {6, 7, . . . , N intI } . (11) 

The remark here is that recursive spin update of intranet is performed asyn- 
chronously. Two-neuron output of intranet has been considered: k G {1,2}, [see 
Eq.©]. The agent's sell buy order 

S(i)^l[s(i,l) + s(i,2)} (12) 

can be identified by other agents. The factor 1/2 normalizes the state space 
{ — 1, 0, 1} of S(i). The 5 = state is interpreted as a passive^^. The impact of local 
Snn(i) and global m in Eq. fTO)) may be interpreted as Keynes' beauty contest US 
according which the stock market notably reflects the mass psychology. 
The agent's decision sell or buy is given by the procedure: 

procedure U S s(i) 
loop III over TVintr repetitions 

1. pick i 6 T randomly; 

2. calculate S nn (i) according Eq.®; 

3. make settings Eq. (fP0")) : 

4. choose randomly intranet spin k S ri ntr ; 

5. calculate the local field at k according Eq. (fTT|) : 

the stochastic field h s t oc h(i, k) is calculated for each node separately; 

6. update S(i) according Eq.(fT2]); 
end of loop III 

The output of the procedure U ss is S(i) variable which can be identified by other 
linked agents [see Eq.®]. 
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3.2. Fitness F(i) concept 

Biological species interact with each other in ways which either increase or decrease 
their fitness. The concept is generalizable to survival of strategies undergoing se- 
lection of the economic agents. The extremal dynamics [see Eq.Q] governing co- 
evolution is based on the knowledge of the local fitness F(i) that expresses an abil- 
ity/inability to survive in the competitive environment. Th e selection and survival 
according fitne ss has be en used for spin market models The fitness-dependent 
link formation I 22 | 24 | 25 | j g conce ptually close to our present formulation. In our 
model, the local fitness is defined as the integral over the history of agent's gains 
and losses calculated from market situation and sell buy orders 

F( t+1 '(j) = + [-c roW + ^NWfO, 1)] . (13) 

The formulation is based on the combination of minority game profit and external 
influences. The relationship of individual spin S and majority is qu antified by the 
— CoS^\i)m^ term with cq > 0. The impact of exogenous news I^SI is included via 
random Gaussian term (white noise) c ran d-^*' ) (0, 1) that is common for all agents. 



3.3. Adaptivity procedure UAd(ia), herding 

Many mathematical models in the social sciences assume that humans can be de- 
scribed as "rational" entities. Nevertheless, the most people are only partly rational 
and in fact are emotional/irrational in t heir decisions. This property is incorporated 
into the concepts of bounded rationality^^ and herding of agents. As a concrete ex- 
ample can serve the followers with their believe that imitation of given strategy 
owned by the social neighbors (selected with care according to their fitness) would 
bring a future benefit to them. The assumes the transfer of information in format 
IIj along the edge of extranet. The adaption of i a starts with the random pick of 
n& G /out that checks a prototype node i pro t = A„ a (i a ) among the nearest neighbors. 
It should be stressed that links are placed by requesting fitness preferences, there- 
fore, the adaptation is random only in the sense that it is not directed to specific 
n a . 

The adaption of IIj(i a ) to nj(i prot ) is described by updates 

Jintr(ia, k, q) <- Jintr(«a, fc, - r/) + ^intr («prot , fc, q) , (14) 

where plasticity parameter r/ £ (0,1) expresses how quickly the follower i a learns 
a strategy. The repeated application of f/Ad(i a ) yields entropy consumption and 
strategic uniformization of market that works against distinctive thinking. 



3.4. Reconnection rules U^i e (i r ), network dynamics 

Recently, the interest in complex networks has been extended to the seeking for 
the local rules governing the build-up of social and technological networks. Several 
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principles have been exploited for this purpose. As an example we mention t he net- 
work that shows marks of age^^Hlor inter- agent communication across the net 1 ^ 1 30 ] 
The core of these methodologies rely on the particular mechanism of preferential 
attachment suggested by Barabasi and Albert 1 ^. A principal distinction from men- 
tioned work is that our present study directs attention on stationarity conditions 
and constant vertex number. 

Let us turn attention to our former proposal of centralized single leader model^^. 
Intuitively, it is implausible that L — l followers can continuously identify the leader 
and its strategy can be adapted by them in the limit of very large L. The reason 
for the revision of single leader picture is that the larger the market is, a more 
demanding and time consuming is the technical analysis of the follower that claims 
to localize socially attainable local leader. Accordingly, we are interested in a more 
general multi-leader stationary society influenced by a varying group of irregularly 
distributed local leaders. 

This paper aims to suggest indirect method of generation of segmented market 
which includes many competing local leaders. The network is reconstructe d indi- 
rectly via the knowledge extracted from the random walk process on the we ^l 32 | 33 ] 
The walk is carried out by the " assistant/informant" agent several times emitted 
from the same source node i T . The informant behaves like the Web surfer read- 
ing pages, jumping from one to another by clicking randomly on web links. To 
be more precise, we define the notion of repeated random walk as a set of linked 
nodes RRW(i r , N Icp , Ad op th) C L. The set is obtained by performing Af pat h steps 
that are N rcp times repeated from the origin i T occupied by the agent seeking for 
a nearby local leader. The fitness ranked by walker - informant allows to propose 
preferential attachments that are expected to yield formation of segmented market 
with impediments to the free flow of information. The procedure of reconnection to 
a local leader denoted as Un. e (i T ) is build upon the sub-procedure of edge pruning 
£^Re w (*n*B) that calls the best connection proposal U^ eB (i T ). The sub-procedures 
are specified in bellow: 

procedure C/R eB (*r) = *b 

(1). loop IV over N rcp repetitions 

that start from the initial condition i\ = i T 
loop V the execution of Adepth iterations 

k+i — X ni {ii) , I = 1, . . . , A/depth , (15) 



for random choices links ni G Tout 
end of loop V ; 
end of loop IV ; 

(2) . the comparison of A^opth x A^ rep values of F(ii) 

(3) . and localization of the agent zb according 



F(i B ) = . max Fin) (16) 

iiERRW 
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t G£ Ai V(Ai)V» 

(a) (b) 

Fig. 1. (a) The time evolution of the log-price returns defined as magnetization mS l \ 
Part (b) shows the scaling properties of the pdf of cumulated returns = nS t+l ^ + 

m (t+2) _|_____|_ m (t+At) p arame ter a — 1.3 determined at first from the fit well coincides 
with the scaling collapse of eight distributions obtained for trading time separations At = 
1,2,. ..,8.). 

via the set RRW(j r , AT rep ; iVdcpth) bounded by the radius iVd ep th 
formed by nodes ii, *2) - - - > iN d<vth visited according 
loops IV, V. 

The output of C/R CB (i r ) is agent «b £ T that is the candidate for future connection 
from node i T . 



The pruning is performed according: 

procedure L/r Cw (^ ^b) = Ur Cv (h-, Ur Cb (»r)) = "w 

loop VI over the N out nearest neighbors. Determine the "weakest" 
(worst) connection n w with the smallest fitness 

F(X n Ji T ))= min F(X n (i T )) (17) 

nel out 

within the actual existing connections 
end of loop VI . 
The output of U^ CB (i r ) is the index n w G I ou t- 

Finally, if F(X n ^(i T )) < F(«b) the output of E/R e (*r) is the update 

*n w (ir)<-tB (18) 

conditioned by the requirement that no multiple connections between i r and are 
established. It also forbids self-connection loops (X n (i T ) = 

3.5. Extremal dynamics t/Ex(*minF)> entry and exit of market 
The extremal dynamics l ^4 | 35 | 36 ] wnerl interpreted as a bankruptcy of firm or exit 
of strategy from market could be considered to be the principal mechanism of the 
economic co-evolution. The idea of Bak-Sneppen model can be easily converted to 
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1000 




(b) 



Fig. 2. (a) The autocorrelation function of the logarithmic price return mM\ (b) The 
indication of the power-law long-time regime of the autocorrelation function of log-price 
volatility |m (t) |. 




(b) 



Fig. 3. (a) The power-law pdf(fc( m )) with local effective exponents 7^ n ' = —1.07 (when 
fc< in ) > 20), and 7 ( in ) = -1.07 [for k (il ^ € (5,20)]. (b) The Zipf's 7F = 0.3 law of fitness. 




0.075 



:o.o5 



0.025 



(b) 




16 32 64 



Fig. 4. (a) The pdf of clustering coefficients. In part (b) the pdf of the minimum path 
ways between randomly selected pair of nodes is plotted. 



strategic variables. In that case the least fit strategy developed by z m i n F is replaced 
by a random candidate from certain limited strategic space. No reconnections are 
assumed within the procedure C/ex itself. This simplification is justified by the as- 
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(a) 




Fig. 5. (a) The snapshot of adjacency L x L matrix (the point denotes the connection). 
The cyclic L-gon maps onto the matrix diagonal, (b) The subgraph {i 6 r,fc' m -'-' > 50} 
(the larger circle belongs to larger k^ ln \ (c) The epochs of topology monitored by means 
of selection of highly preferred {i, i € F, fc' m ) (i) > 100} C I\ The members of such group 
may be vaguely identified as local leaders. 

sumption that low fitness agents are rarely attached [i.e. A;( m )(i m i n F) is relatively 
small] . This reluctance against attachment follows from the preferences incorporated 
into [/Re (*r)- The extremal event in the present formulation means that instant value 
of the strategic variable is immediately replaced by Gaussian distributed random 
number N(0, f{...}) of dispersion &{...}• The new strategy enters thought 

Jintr(«minF, k, q) «- N(0, (Tj intI ) , (19) 

^stoch(iminF, k, q) «- N(0, a hstoch ) , 
F(i minF ) ir- N(0,a F ). 

These updates necessitate the fixing of three independent dispersion parameters 

4. Simulation results 

4.1. Selection of parameters 

The statistics as presented in the next subsection is determined by the choice of 13 
free parameters. Their adjustment is a nontrivial task as whenever in the modeling 
of agency. To attain at least qualitative agreement with current statistical concepts 
in finance values (see Table [T]) have been suggested by the optimization 

in the parametric space. A particular requirement is to keep the spin dynamics 
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Table 1. List of 13 numerical parameters for which the statistics is 
presented. 



number of 


extranet nodes 


L = 500 




node outputs 


N out = 10 




random search steps 


-^depth = 6 




repeated net searches 


JVrcp = 6 




intranet nodes 


AW = 8 


dispersion of 


^intr 


^intr = 1 




^stoch 
F 


crp = 0.1 


probability of 


reconnection 


Pr c = 0.01 




adaptive move 


P Ad = 0.2 


adaptivity 


parameter 


n = 0.1 


fitness 


minority game 


co = 1 




news 


Cran = except Tab[2] 



Table 2. The study of unexpected information about eco- 
nomic performance or political situations incorporated into 
fitness. The comparison of averages corresponding to two 
combinations of parameters cq, c ran( j of local fitness Eg. 1131 . 
The measure Pm,4 denotes the probability that tyv^' does 
not alter sign during the four subsequent steps m' 1 ', rn^ t+1 \ 
m (t+2) ^ m (t+3) (f our ti mes bearish or four times bullish 
stock). The table shows that risk aversion accompanied by 
the probable passive states with reduced (|rm')|) can is- 
sue from the high exogenous influence. The calculation of 
PmA indicates that news can break bearish (bullish) se- 
quences. The population of local leaders quantified by ex- 
ceptionally connected nodes characterized by probability 

-Pfe>100 = Efc(in) >100 Pdf(fc( m )). 



CO C rand 


(|m (t, |> Pm,4 


-Pfc>100 (C) (imia) 


1 0.00 
1 0.05 
1 0.10 


0.0138 0.226 
0.0082 0.165 
0.0029 0.091 


0.02 0.55 15.2 
0.02 0.56 18.4 
0.02 0.57 15.2 



much faster than . An additional constraint aims at at tain ing; vicinity of a critical 
regime, where the power-law distributions are obeyed that are much studied 
topics inside Econophysics. Various mechanisms of the generation of power-laws 
have been recently summarized in review article^]. 

4.2. Distributions and averages 

Starting from random initial condition, we let the system evolve at least 5.10 4 Monte 
Carlo steps per node. The data for calculation of averages of interest have been 
collected from the next ~ 10 7 Monte Carlo steps. Their validity is verified for several 
independent trials. The comparison of selected averages corresponding to different 
fitness is given by TabJ21 The table indicates that external news yield more careful 
behavior of agents and preference of passive states 5 = that anomalously sharpen 
the central part of probability density function (pdf) of m denoted as pdf(m). For 
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given example, the impact of c rall on the topology is marginal. 

In further we concentrate on the case c ran d = 0. The statistical treatment (see 
FigfTJ) leads to the fat tailed pdf(m) that have been fitted in particular by the power 
law asymptotics of the realistic effective exponent a ~ 1 .3 12JJ The dynam- 

ical features of the price dynamics are highlighted by the autocorrelation functions 
in FigfSJ The clustering of the volatility \m^ | of the log-price returns is observed. In 
that case the power-law indicates the occurrence of the long time memory. Despite 
of the fact that exponent —0.64 does not strictly reproduce empirical findings —0.34 

and — 0.2^311, the recovery of log-slope from (—1,0) range is quite encouraging 
for the perspectives of the model. The situation is even more interesting since the 
power-law pdf(F) [see Fig|3ljb)] histogramed for F > can be interpreted as Zipf 's 
7f = 0.3 law for survival abilities of strategies. 

Let us turn attention to the issues of network statistics depicted in FigsEl^a) J4] 
and0 The node degree k^ ln '(j) = J^ier ^2nei t $j,x n (i) accounts for incoming links 
of node j. The stationary regime generates sequence of networks with broad-scale 
pdf of node degrees. By fitting of pdf(k^) we have identified two partial effective 
exponents 7 in ~ -1.07 [for fc( in > 6 (5,20)], and -1.72 (as fc( in ) < 20) from the 
law pdf(fc m ) ~ [fc( m )] 7 . Several real exponents are provided here for illustrative 
purposes. The value 1.81 ^3 corresponds to the collection of e-mail addresses. The 
exponent —1.2 belongs to the coauthorship network SB. 

I n fu rther the organization of the network is done using clustering coeffic ient 
C(i)SHthat is a local measure of interrelatedness of triplets or social transitivity^^. 
For directed network the tendency can be measured by C{i) — e(i) / (N out (N out — 1)), 
where e(i) = E^^^i"" "' S x„ 1 {Xn 2 (i)),x„ 3 ({) stands for the number of the links 
between neighbors X ni (X n2 (i)) and X n3 (i) attained from i. The maximum number 
of links N out (N out — 1) normalizes the expression for C(i). The object of interest is 
the mean (C) . As usual, it is meaningful to compare the mean clustering coefficients 
of two distinct network reconnection modes. For the network with randomized / ou t 
links we obtained (C rand ) ~ 0.02, while (C) ~ 0.55 ((C*)/(C cand ) ^ 27.5) reached 
by C/ro is much higher and thus empirically more relevant . We have computed 
the average of the minimum path way for the partially random net (^min.rand) — 
2.9 (with the circular subgraph untouched). The action of U^, e (i T ) provides (l m i n ) — 
15.4 however much smaller maximum iniin,max = 3 of pdf(Z m j n ). The combination of 
above topological attributes supports both small-world and scale-free concepts of 
network statistics. Fig|5] indicates that RRW(., ., .) invokes generation of modular- 
topology of several local densely connected leaders. 



5. Conclusions 

In this paper we reveal some interesting properties of the model of spin dynamics 
on a complex evolving network. We demonstrated that our model is relevant for 
the description of stock market statistics. The stationary regime is formed due to 
balance between information entropy inflow produced by the extremal dynamics 
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and stochastic outer sources compensated by the entropy outflow caused by the 
adaptive moves. 

We are conscious that minority game term —CQS^(i)m^ of fitness only roughly 
describes the financial profits at th e re al stock markets. In further research we plan 
the comparison with $-game model S3 that represents a step closer to stock market 
reality. Additional improvements of fitness could be done by considering the impact 
of regulatory and non regulatory halts and delays. 

The most emergent aspects of agent systems near the criticality are the power- 
law distributions of the topological and spin characteristics that occcur as a conse- 
quence of self-organization processes. The necessary remark in this context is that 
exponents of such dependencies are non-universal, i.e. they can vary from one set 
of free parameters to another. 

The adjustment of parameters can be formulated as multi-objective optimization 
computationally demanding task. Due to its comprehensive features, the problem 
is planned to be discussed in future works. 

The site http://158.197.33.91/%7Ehorvath/selfstructured_net/ provides our 
simulation C++ code and supplementary materials. 
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